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Abstract 

A localized charged particle oscillating near a reflecting boundary is considered as a model for 
non-cancellation of vacuum fluctuations. Although the mean velocity of the particle is sinusoidal, 
the velocity variance produced by vacuum fluctuations can either grow or decrease linearly in time, 
depending upon the product of the oscillation frequency and the distance to the boundary. This 
amounts to heating or cooling, arising from non-cancellation of electric field fluctuations, which 
are otherwise anticorrelated in time. Similar non-cancellations arise in quantum field effects in 
time-dependent curved spacetimes. We give some estimates of the magnitude of the effect, and 
discuss its potential observability. We also compare the effects of vacuum fluctuations with the 
shot noise due to emission of a finite number of photons. We find that the two effects can be 
comparable in magnitude, but have distinct characteristics, and hence could be distinguished in 
an experiment. 
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I. INTRODUCTION 



Consider a localized charged particle coupled to quantum electromagnetic field fluctu- 
ations in the vacuum state. We will treat it as a classical particle, but more generally it 
can be viewed as a quantum particle in a wavepacket state sharply peaked in space. Be- 
cause the vacuum is the state of lowest energy of the quantum field, the particle cannot, 
on average, acquire energy from the electromagnetic field. This does not prevent energy 
fluctuations which are within the limits set by the energy-time uncertainty principle. The 
particle can acquire additional energy from an electric field fluctuation, but the energy must 
be surrendered on a timescale inversely proportional to the magnitude of the energy. Energy 
conservation is enforced by temporally anticorrelated electric field fluctuations, which are 
guaranteed to take back the energy within the allowed time. Thus on the average, neither 
the particle nor the quantum field gains energy. 

This holds in any static situation, including one where reflecting boundaries are present. 
Although classical image charge effects can be present, no net energy may be extracted 
from the vacuum. A model with a charge maintained at fixed mean distance from a plane 
mirror was treated in Ref. pp. Switching on the effect of the mirror can cause the particle's 
mean squared velocity to either increase or decrease, but after transients have died away, it 
approaches a constant. This need not be the case in a time-dependent situation, which will 
be the topic of this paper. The cause of the time-dependence may be a source of energy, 
so it is now possible for the particle's energy to either grow or decrease in time. However, 
one may also view the time-dependence as upsetting the anticorrelated fluctuations which 
are present in a static situation. In the static case, the anticorrelated fluctuation takes 
exactly the amount of energy obtained by the particle in a previous fluctuation. The time- 
dependence may either enhance or suppress the magnitude of the the second fluctuation, 
resulting in either a decrease or increase, respectively, of the particle's energy. We will see 
both possibilities illustrated in the model discussed in Sect. [TTj 

Examples of non- cancellation of field fluctuations arise in cosmology. One is Brownian 
motion of charged particles in an expanding universe [2]. Other examples were discussed 
in Refs. [3H5], where it was argued that quantum stress tensor fluctuations during inflation 
can lead to density and gravity wave perturbations which depend upon the total expansion 
during inflation. In the present paper, we consider a simple flat space model which is of 
interest both in its own right, and as an analog model for effects in curved spacetime. 
Lorentz-Heaviside units with c = h = 1 will be used. 



II. THE MODEL 



A. Formulation and Calculations 

Our model consists of a particle of mass m and electric charge q undergoing bounded, 
non-relativistic motion in a direction normal to a perfectly reflecting plane mirror. We take 
this to be the z-direction, and write 

z(t) = d + Af(t), (1) 

where d is the mean distance to the mirror, A > is the amplitude of the motion, and f(t) 
is a dimensionless function which we later take to be sinusoidal. We require z(t) > for all 
t and |i(t)| = A \f(t) \ <C 1 We assume that the components of the particle's velocity satisfy 
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a Langevin equation, 

Vi = ^- Ei(x,t) , (2) 
m 

where x = x(£) is the spatial location of the particle at time t. Here E is the total electric 
field, including both a classical applied field, including possible image charge effects, and the 
quantized electric field. This is the usual equation of motion for a non-relativistic charged 
particle when magnetic forces are neglected. Our key assumption is that it may be used in 
the presence of a fluctuating electric field. For now we ignore dissipation effects, which have 



been discussed in Refs. [HI U\. We will treat dissipation by emitted radiation in Sect. Ill A 
Note that an alternative to moving the charge with the mirror fixed is to move the mirror, or 
to use a charge moving at constant speed near a corrugated mirror. The latter strategy was 
first used by Smith and Purcell [H] to create radiation, and is the basis of the free electron 
laser. 

With the initial condition Vi(t ) = 0, we may integrate the Langevin equation and then 
take expectation values in the electromagnetic field vacuum state to write the variance in Vi 
as a double time integral of the electric field correlation function: 

(A^ 2 (t)> = 1- [ & x f dt 2 [(£?i(xi,ti) Ei(x 2 ,t 2 )) - (Eiixuh)) (Ei(x 2 ,h))} . (3) 
m J t J t 

Here Xi = x(£i) and x 2 = x(t 2 ), the spatial locations of the particle at times t% and t 2 , 
respectively. Any classical part to the electric field will cancel in the correlation function. 
For now, we focus on the quantum part of the electric field, for which (^(x, t)) =0. We are 
interested only in the effect of the boundary, as the empty space correlation function will 
not produce any growing terms in (Avf(t)). The quantum electric field correlation function 
may be written as a sum of an empty space part and a boundary correction. We drop the 
former and write 

2 ft ft 

(Av*) = ^ dt x \ dhiEifatjEifah^b, (4) 

m Jto Jto 

where the subscript b indicates the boundary correction to the two-point function. These 
corrections may be found by the method of images, and are [H] 

(E x (^, tl )E x ^ 2 ,t 2 )) b = - 2 [o + 7^;\ 21 3 (5) 
for a transverse direction, and 

{E z {*lM)E z {* 2 ,t 2 ))h = 7r 2[ r 2_(^ +;Z2 )2]2 ( 6 ) 

for the longitudinal direction, where r = t% — 1 2 , and z\ = z(t±), ect. Here we assume that 
the particle does not move far compared to the distance to the mirror, and have equated 
the coordinates in the transverse directions, x\ = x 2 and y\ = y 2 . Note that, for example, 
(v x v z )b = 0, so there will be no correlation between the random motion in the transverse 
and longitudinal directions. 

Next we assume that d, and Taylor expand the two-point functions to second 

order in A. The integrand for the longitudinal variance becomes 



T 
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-{2d + A[f{t l ) + f{t 2 Wf 

1 + 7 5 8d L, ^/(*i) + m)] + 2 K ± ^ Am) + f{t 2 )f . (7) 



(r 2 - Ad 2 ) 2 (r 2 - Ad 2 ) 3 LJ K 1J J K in (r 2 - U 2 \ 
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We are seeking contributions to (Au?) which grow in time. The zeroth order term describes 
the case of a stationary charge, which was treated in Ref. [T|, and gives a constant contribu- 
tion. The first order term yields a purely oscillatory function when f(t) is sinusoidal. Thus 
we omit both of these terms and focus on the second order term. Note that the r-dependent 
part of this term may be written as a total derivative 



2(r 2 + 20d 2 ) d 4 d 2 d 2 



T 



r 



The function G(r) may be expressed in terms of logarithmic functions, but we will not need 
its explicit form, beyond the fact that it has only a logarithmic singularity at r = 0. 

Now we assume that f(t) and its first three derivatives vanish in the past and future. 
This allows us to integrate over all t\ and ti, and to perform integrations by parts with no 
boundary terms. Thus we may write 

/oo poo o2 a2 

dt 1 dt 2 F z (r)[f 2 (t 1 ) + f 2 (t 2 )} = J dt 1 dt 2 G(r)-^^[f 2 (t 1 ) + f 2 (t 2 )]=0. (9) 

This implies that only the cross term in the last term in Eq. (|7|) can give a nonzero contri- 
bution. Now we may write 

9 n 2 rco poo 

(Av 2 z ) = -^A 2 / / dh dt 2 F z (r) fit,) f(t 2 ) . (10) 

7T ™ J-oo J-oo 

Next we adopt a specific form for f{t\), which is f{t\) = sin(a;£i) for < tx < t and 
f(ti) = for t\ < and t\ > t. The approximate signs indicate that / should fall smoothly 
to zero at the end points of the interval. This describes a charge which oscillates sinusoidally 
at angular frequency u for a time t. This sinusoidal motion could be driven by a classical 
electric field of the form Ef(t) = —E sin(atf), in which case 

A= q A- (ID 
muj z 



The integration in Eq. (10) is effectively over a square of side t. Next, we change integration 
variables to r and u = t\ + 1 2 . Because F z (r) falls to zero rapidly if |t| 3> d, and because 
we assume t 3> d, the integration on r may be taken over an infinite range. However, the u 
integration is restricted to a finite interval: 

i a 2 r 2t r°° 

(Av 2 z ) = ^~ 2 — 2 A 2 \ du / dTF z (r)[cos(cor) - cos(um)] . (12) 

Z7T m Jq J^oa 

The integral of the cos(uu) term will generate an entirely oscillatory contribution, which 
may be ignored compared to the linearly growing term, so we may write 



2 a 2 

(A^> w 4 —2 A " t 
n 2 m 2 



, , , 00 (r 2 + 20d 2 



{T 2_ Ad 2 Y ^(ur) dr 



(13) 



At this point, it is useful to note that r should have a small, negative imaginary part in 
Eqs. Q and (§)• This arises because these two-point functions are expressible as integrals 
of the form 



oo 

3 „— iuJT 



duu d e- lulT , (14) 



o 
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which are absolutely convergent if Im(r) < 0. We can implement this condition by replacing 
r by r — it in Eq. ( 13 ), where e is a small positive real number. We can write the denominator 
in the integrand as 



[(r - it) 2 - 4rf 2 ] 4 = (r - it + 2d)\r -it- 2d) 4 , (15) 

revealing that there are two fourth-order poles in the upper half-plane at r = ±2d + ie. Next 
we write cos(o;t) in terms of complex exponentials. The r integration is along the real axis, 
so the e~ luJT term gives no contribution when the contour is closed in the lower half-plane. 
The e tU)T term yields the residues of the two poles when the contour is closed in the upper 
half-plane. The sum of the residues is a real function. 



B. Key Results 



The result of the evaluations of the longitudinal velocity variance, after using Eq. (11 ), is 

^,4 7712 

167rm 4 d 



(^) = 43l^t, (16) 



where 

R. = 

2£ 



R* = ^[(3 - 5£ 2 ) sin(20 + 2£ (£ 2 - 3) cos(2£)] , (17) 



and £ = uj d. 

The same mathematical technique holds for the transverse direction; only the precise 
form of the integrand changes. Let F z — > F x , where 

4(40rf 4 + 34rfV + r 4 ) 
Fx{t) ~ [(r-ie)*-4d2]5 ' (18) 

In this case, there are two fifth-order poles in the upper half-plane, but otherwise the evalu- 
ation procedure is the same. Now the velocity variance in the x-direction, which is also the 
mean squared velocity in this direction, is found to be 

< A "*> = <»*> = 10^- < 19 > 

where 

= ^FT^[( 4 £ 2 - 3 ) sin ( 2 ^) + 6 £ co < 2 0] , (20) 

Note that the Ri, which are dimensionless, are proportional to the rate of change of the 
corresponding velocity variance: 

16nm 4 d d(Avf) 

m) = ^ET^T- (21) 

These quantities are illustrated in Fig. [Tj 

Of significant interest here is that for both the longitudinal and transverse components, 
the coefficient of the time dependence of (Av f) can be either positive or negative, depend- 
ing on the frequency of the oscillation and distance to the mirror. These results can be 
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FIG. 1: This graph shows the dimensionless quantities R z and R x , which are proportional to 
the rates of change of velocity variance in the longitudinal direction and a transverse direction, 
respectively. Here £ = u d. 



interpreted in terms of non-cancellation of previously anticorrelated electric field fluctua- 
tions. When there is linear growth, the fluctuations are adding energy to the particle on 
average. Similarly, a linear decrease signifies that they are removing energy, which could be 
described as a "cooling mode" . The latter effect can only go so far, and at some point our 
approximation of localized particles would break down. 

It is also of interest to examine the low and high frequency limits of the above results. 
At low frequency, £ <C 1, we have 

and at high frequency, £ ^> 1, 

^ ~ Sin ^ ' < Aw 2> ~ #4^7 c°s(2£) . (23) 

lo7rm 4 a 167rm 4 a£ 

Note that the effect tends to be larger in a transverse direction than in the longitudinal 
direction, especially at high frequencies. 

Next we wish to make some estimates of the magnitude of the heating or cooling effect. 
We do this by defining a change in effective temperature for the i-direction, ATj, by 

^m(Av*) = h B AT l , (24) 

where k B is Boltzmann's constant. Strictly speaking, this is not a real temperature, since it 



is not isotropic, but it is a useful measure of the size of the effect. From either of Eqs. (16) 



or (19), we find 

1 ^0 



1 167ik B m 3 d 1 ^ ^ 
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This may be expressed as 



— 10 - 8 -(iwb)(^)fc) fl " (26) 

where we have replaced Eq/2 by /, the power per unit area in a plane electromagnetic wave 
with peak electric field E . We have also set q = e, the electronic charge. 

Our approximation of a perfectly reflecting plate should hold both for modes whose 
wavelength is of order d and at angular frequencies of order u. Note that £ = 2ird/\, where 
A is the wavelength of the driving field. From Fig. [TJ we see that R z reaches its maximum 
value of about 0.5 at £ ~ 2.5 and R x first reaches its maximum of about 1.0 at £ ~ 4. Both 
of these correspond to A > d. If d > Ap, the plasma wavelength of the metal in the plate 
which can be in the range of 0.1/im, then the perfect reflectivity assumption should be valid. 
Ultimately, whether this effect can be measured in a realistic experiment depends upon the 
sensitivity of temperature measurements, the power intensity / of the driving field which can 
be used, and the time t which can be achieved. On the latter point, it is possible that planar 
Penning traps will be able to achieve very long coherence times with single electrons |10j . 

As noted earlier, (Av^) = (v%) because the mean transverse velocity vanishes, (v x ) = 0. 
Thus the increased drift in the transverse directions when (v%.) > is a signature of this 
effect. When (v%.) < due to the shift in electromagnetic vacuum fluctuations, then we need 
to interpret the effect as a reduction in mean squared transverse velocity, with a positive 
contribution coming from other effects, such as quantum uncertainty in speed, classical 



thermal effects, or shot noise (to be discussed in Sect. Ill B ). This reduction is closely 
related to the phenomenon of negative energy density in quantum field theory, whereby it 
is possible to reduce to local energy density below the vacuum level with either boundaries 
or quantum coherence effects 

In the longitudinal direction, there is a nonzero mean velocity given by the response to 
the classical driving force. The time averaged square of this velocity is 



<*>'=2 =2^» ' < 27 > 



It is of interest to compare this quantity with the quantum variance given by Eq. (16), and 
write 

Given that we expect £ ~ R z ~ 0(1), it is possible for the random motion produced by 
electric field fluctuations to approach a significant fraction of the classical motion. 

So far, we have treated the charges as classical point particles, but the same conclusions 
about changes in (At> 4 2 ) should hold for quantum particles if they are localized in space on 
scales small compared to d, the distance to the mirror. Ehrenfest's theorem tells us that 
classical equations of motion hold, on average, in quantum mechanics. We are concerned 
here with changes in the variance of the particles' momentum, which follow from momentum 
conservation considerations. 
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III. RADIATION AND SHOT NOISE 



A. Radiated Power 



In this subsection, we will address the dissipation effect due to emitted electromagnetic 
radiation. First, we determine the average power radiated by our oscillating charge. It acts 
as an electric dipole, and so we may use the far-zone field formulas, with the method of 
images to obtain the field components. (This problem has been extensively studied in the 
literature. See, for example, Ref. [12] for further references and a detailed treatment of the 
near-zone fields.) What follows are the non-zero field components at a point of observation 
P located in the x — z plane, where r\ is the distance from P to the real dipole, and r 2 the 
distance from P to the image dipole. We have assumed that P is far enough away that both 
r% and r 2 have approximately the same polar angle 9. First, the z-component: 

sin 2 fl , fe iuri e ur2 \ . . 

E z = — Pe u 2 — + — ) (29) 



4/T \ n r 2 

where p e is the peak value of the oscillating electric dipole moment and u is the frequency. 
From here, we make further approximations: given a distance 2d separating the dipoles, we 
can let r\ ~ r + dcos6 and r 2 ~ r — dcosd. Further, since we are assuming d <C r, we 
approximate r\ ~ r 2 w r in the denominators. The z-component is then 

p e co 2 e iu)T 

E z = sin 2 9 cos(uid cos 9) (30) 

2n r 

Similarly, for the other non-zero field components, we have: 

p e u 2 e iwr 

E x = — sin 9 cos 9 cos(ud cos 9) , (31) 

Air r 

and 

p e to 2 e tujr 

Hy = — ^ sin^cos^c/ cos9) (32) 

The next step is to obtain P(9), the power radiated per unit solid angle in the direction 
of a unit vector n = sin^x + cos#z. From the Poynting vector, we find 

P(9) = r 2 n • (E x H*) = r 2 (sin^x + cos# z) • (—E Z H* x + E X H* z) 
p 2 UJ^ 

[sin 9 cos 2 (cud cos 9) + sin 2 9 cos 2 9 cos 2 (cud cos 9)} 



8tt 2 

p\ ^ 



sin ^cos (ujd cos 9)). (33) 



8tt 2 

We next integrate P(9) to obtain the total power radiated: 



t*2n / > 7i"/2 

P T = I / P(9)dn. (34) 



'0 

Let u = cos 9 and use p e = qA and £ = co d to write 



Pt = ^ /V - u 2 ) cos^u)du = ^{s+ £[-2*cos(2£) + sm(20]} = , 

(35) 
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where 

3 



S t = 1 + [-2C cos(20 + sin(2£)] . (36) 



This gives us Pt, the energy radiated per unit time. 

We can write the energy radiated per oscillation cycle, E c , as 

2irP t 1 2 2 , 



E ° = — = *<l A " S T . (37) 
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The ratio of this quantity to the particle's average kinetic energy is 

E c q 2 A 2 u 3 S T 2q 2 A 2 u 3 S T 



(38) 

St ■ (39) 



(Eian) 3m(v 2 ) 2>mA 2 uj 2 

E c _ 2q 2 u 
(E kiQ ) 3m 



The function St is of order one when £ is of order one. Then, inserting the charge and mass 
values for an electron, as well as our typical frequency value of 10 l4 Hz, the estimate comes 
out to 

jp- « 8 x 10- 9 (40) 

Thus, the electron radiates only a few parts per billion of its own kinetic energy per 
cycle. The small value of this ratio shows that the electron with our driving field is a weakly 
damped driven oscillator that needs only minimal energy restoration for preservation. The 
emitted radiation is the primary irreducible source of dissipation. This estimate indicates 
that it is reasonable to neglect its dissipative effects on the motion of the paricle. 



B. Shot Noise from Photon Emission 



However, there is another effect arising from the emitted radiation to be considered. 
Because the power radiated by the particle consists of discrete photons, there will be a 
statistical uncertainty in the momentum lost by the particle. This will lead to an additional 
contribution to (Av 2 ), the velocity variance of the particle. Any experiment which seeks to 
measure the effects of vacuum fluctuations on the variance, Eqs. (16) and (19), will have to 



contend with this shot noise as a background. Let Pi be the average power radiated by the 
particle in direction i. Then in time t, an energy and magnitude of momentum of pi = Pit 
will be radiated in this direction, corresponding to a mean number of photons of N{ = Pi t/oj. 
The statistical uncertainty in this number is y/Nl, assuming that the emission of different 
photons are uncorrelated events. This leads to an uncertainty in the i-component of the 
particle's momentum of order 



Api = u \/Wi = ^/Piiut . 
and a variance in the velocity in direction % of 

Put 



Avi 



(41) 



(42) 



where the "s"-subscript refers to shot noise. 
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Now we find the total power radiated in the z-direction. This quantity is found by 
projecting onto the z-axis, and integrating over a hemisphere: 

p2n /*tt/2 /"^/S 

P z = r 2 I I P(9)cos6dtt = 2nr 2 / P(6) cos6 d(cos9) 



JO 
2, ,4 /•! 



u(l - u 2 ) cos 2 { q u)du, (43) 



47T 

where u = cos 9, as before. The result is 

Pz = oif [ " 3 " ^ + {3 ~ C ° S(2 ° + m " + 3sin ( 2 ^)] ' ( 44 ) 
Next, introduce substitutions for the dipole moment as follows: 

p 2 e = q 2 A 2 = q 2 (^) 2 . (45) 

We now have 

P = qAE o ji UQ) 

where, 

j = ^[-3 - 2£ 2 + (3 - 4£ 2 ) cos(20 + 2£(£ 3 + 3 sin(20)] . (47) 
Consequently the mean square velocity in the ^-direction from shot noise is 

Av 2 z = 1 . (48) 



Now compare this effect to that of the electric field fluctuations, using Eq. (16) to write 

(At;*) 4R Z 



Avl S x 



(49) 



Figure [2] compares these effects, showing the relative magnitudes of what would be seen 
without and with quantum electric field fluctuations, as a function of £. 

We can make a similar calculation for the power radiated in the x-direction, and find 

p2tt /•it/2 o 4 rp2 q 

P.^l I PMsi.^A^f. (50) 

Here 

S. = £ [2^® + l] , (51) 
and J2 is a Bessel function of the first kind. We find the x-direction velocity variance to be 
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FIG. 2: The relative magnitudes of velocity variance in the z-direction with only shot noise, S z , 
and with both shot noise and quantum electric field fluctuations, S z + 4R Z . 




FIG. 3: The relative magnitudes of velocity variance in the x-direction with only shot noise, S x , 
and with both shot noise and quantum electric field fluctuations. 



The ratio of the effect of electric field fluctuations to that of shot noise for the transverse 
direction is 

Av 2 sx 3nS x ' { ' 

The same graphical comparison as for the z-direction leads to Fig. [3] We see that the effects 
of quantum electric field fluctuations and of shot noise are comparable in order of magnitude 
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when £ is of order one. However, the sum of the two effects always seems to lead to a positive 
velocity variance. In the limit that £ 3> 1, we find 



Av 



and 

(Aj) 8sin(2Q 
Ai£ ~ 3tt£ 



(55) 



Thus, in the limit of high oscillation frequency or large distance to the mirror, the shot noise 
effect dominates. 



IV. SUMMARY 

In summary, we have presented a model in which charges, such as electrons, moving in 
the quantum electromagnetic vacuum near a mirror may increase or decrease their velocity 
variance. The ultimate energy source is the driving field, but the mechanism can be viewed as 
non-cancellation of anticorrelated electric field fluctuations. The effect is a form of squeezing 
of the particle's velocity uncertainty by the electromagnetic vacuum fluctuations. The most 
striking aspect of this effect is that the mean squared velocity can decrease, corresponding 
to an effective cooling of the charges. Although the effect is normally small, it might be 
observable. 

In our model, we have assumed that the charges move and the mirror remains station- 
ary. However, for non-relativistic motion, one would obtain the same result if the opposite 
were true. A rapidly oscillating mirror is more difficult to achieve, although rapid electrical 
switching of the reflectivity of a mirror might be possible, and has been explored in the 
context of the dynamical Casimir effect, the quantum emission of photons by a moving 
mirror [131 E]. This effect seems to have been recently observed in the context of supercon- 
ducting circuits [15]. Although the effect discussed in the present paper involves exchange 
of kinetic energy between charges and a quantum field in the presence of a boundary, rather 
than quantum creation of photons, it can be viewed as a variant of the dynamical Casimir 
effect. In the latter case, the kinetic energy of the boundary is converted into photons. In 
the model of this paper, it is converted into random motion of a charged particle, but both 
are effects in quantum field theory. 

An alternative to switching of a mirror is the use of charges moving near a corrugated 
mirror, as in the Smith-Purcell effect [5] . In this configuration, the effect studied here should 
also arise. 

We compared the effects of electromagnetic vacuum fluctuations with shot noise due 
to emission of a finite number of photons. The two effects can be of the same order of 
magnitude, but have distinct signatures, so it should be possible to distinguish them exper- 
imentally. 

The effect studied here is also of interest as an analog model for quantum effects in 
cosmology. A curved background spacetime can also cause non-cancellation of otherwise 
anticorrelated fluctuations. Thus the effect discussed here bears some relationship the effects 
studied in Refs. [2H5]. 
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